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ABSTRACT 

Pn-dh’h-d turbine vane heat transfer for rough sur- 
faces was compared with experimental data for both 
vanes and rotors. For the vane comparisons, inlet pres- 
sures varied 1m tween 0.2 and 1 atm., and exit Mach 
mini I** ranged between 0.3 and 0.9. Thus, while a 
sing!* rouji surface vane was used for the tests, the 
effect i\* roughness in wall units varied by more than 
a factor <4 on Comparisons were made for both high 
and low fr* * >tr« am turbulence intensities. For the ro- 
tor. comparisons were made at two Reynolds lumbers 
for each of two inlet flow angles. Results are shown for 
both smooth and rough rotor blades. Two-dimensional 
Navi* r-Stoh' s heat transfer predictions were obtained 
using tin rod* RYCQ3D. Results were obtained using 
both algebraic and k — turbulence models. The al- 
gebras mod* 1 incorporated the Cebeci-Chang rough- 
ness mod* I Hi* /■ — u turbulence model accounts for 
roughness in ih* boundary condition. Roughness causes 
turbulent llo\v over the vane surface. Even after ac- 
counting fnr transition, surface roughness significantly 
increased In at transfer compared to a smooth surface. 
The k — + results agreed better with the data than 
the Cebeci-Chang model. The low Reynolds number 
k — u: model did not accurately account for roughness 
at low freest ream turbulence levels. The high Reynolds 
number version of this model was more suitable at low 
freestream turbulence levels. 


h - Roughness height 

k - Turbulent kinetic energy 

/ - Length scale 

M - Mach number 

Nu - Nusselt number 

Re - Reynolds number 

s - Surface distance 

Tv - Turbulence intensity 

U - Velocity 

y - Normal distance from vane surface 

0 - Relative flow angle 

k - Von Kami an constant 

A - Roughness density parameter 
p - Molecular viscosit y 

Hx - Turbulent eddy viscosity 

p - Density 

uj - Specific dissipation rate 


Subscripts 

EQ 

FS 

i 

IN 

0 

I 


Equivalent 

Freestream 

Inner region 

Gas inlet 

solid boundary 

Value at first grid line 

Vane row exit 


Superscript 

+ - Normalized 


Nomenclature 

- Near wall damping coefficient 
Cj/2 - Friction factor 
C p - Pressure coefficient 
C x - Axial chord 


INTRODUCTION 

The ability to predict the effects of surface rough- 
ness on both turbine aerodynamics and heat transfer is 
important. Several researchers have reported decreases 
in turbine efficiency of up to several points due to 
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surface roughness. Among these are Kind ei al.(199S), 
Boynton et al.( 1 992), and Bammert and Stanst- 
ede( 1972 ,1976). In addition to causing a decrease 
in aerodynamic efficiency, surface roughness can sig- 
nificantly modify turbine blade surface heat trans- 
fer. Dunn et al.(1994), Blair(1994), Tarada and 
Suzuki(1993), and Abuaf et al.(1997) among others 
showed the effects of surface roughness on turbine blade 
heat transfer. The degree to which surface roughness af- 
fected the surface heat transfer varied among the inves- 
tigators. Accurate heat transfer predictions for rough 
surface turbine blades are important in accurately pre- 
dicting turbine blade life. Surface roughness can cause 
a laminar boundary layer to become turbulent, and 
increase the heat transfer coefficient over that for a 
smooth blade. The increased heat transfer coefficient 
may or may not increase the heat load to the blade. If 
the surface roughness is caused by deposition of low con- 
ductivity material, the heat load may decrease. How- 
ever, if the roughness is caused by erosion the heat load 
increases at the same time the blade strength is de- 
creased. The present work is focused on verifying pre- 
dictions for the external heat transfer coefficient, and 
does not address the effect on blade temperature. 

Several investigators have analyzed the heat transfer 
to rough surface blades. One approach, advocated by 
Taylor et al.(1985), and utilized by Tolpadi and Craw- 
ford(1998), and Tarada(1990) is to model the rough- 
ness as geometrical elemen ts attached to a smooth sur- 
face. To account for the blockage of these elements 
the equations of continuity, momentum and energy are 
modified in the near wall region. Another approach 
is to assume the roughness affects the flow in a way 
similar to sand grain roughness. In this approach the 
turbulent eddy viscosity is increased based on corre- 
lations for sand grain roughness. Here the physical 
roughness is related to equivalent sand grain rough- 
ness by empirical correlations. For algebraic turbulence 
models Cebeci and Chang(1978) recommended a pro- 
cedure to modify the turbulent viscosity for sand grain 
roughness. Their model was used by Boyle and Civin- 
skas(1991) to predict heat transfer for rough turbine 
blades. Wilcox(1994) accounted for surface roughness 
by modifying the boundary condition in the k —w tur- 
bulence model. It is the latter approach that is utilized 
in the work reported herein. This approach is simpler 
than the one advocated by Taylor et al.(!985), and it 
was felt that the variation in roughness properties, even 
within a single roughness trace, was large enough to dis- 
courage modeling the roughness as a series of repeated 
geometric elements. 


The work reported herein used the data obtained 
by Boyle et al.(2000) to identify an appropriate means 
of predicting the effects of surface roughness on turbine 
blade heat transfer. It is felt that the approach identi- 
fied would also be a suitable candidate for the verifica- 
tion of turbine aerodynamic performance. In addition, 
comparisons are shown with the midspan rotor blade 
heat transfer presented by Blair(1994). This is done to 
show that the conclusions drawn from the comparisons 
with the data of Boyle et al.(2000), also apply for the 
comparisons with the data of Blair(1994). 

DESCRIPTION of ANALYSIS 

A two-dimensional Navier-Stokes analysis was used 
to predict the effect of surface roughness on vane surface 
heat transfer. The computer code used was the quasi- 
three dimensional analysis described by Chima(1987) 
and by Chima an Yokota(1988). For the test configu- 
ration analyzed there were no three-dimensional effects 
expected. The analysis developed by Chima(1987), and 
Chima and Yokota(1988) is a thin layer Navier-Stokes 
formulation, that achieves a steady-state solution us- 
ing a Runge-Kutta time marching approach. Implicit 
residual smoothing is also used. 

Two algebraic models, Baldwin and Lomax(1978) 
and Chima et al.(1993), and one two-equation k— & tur- 
bulence model, Chima(1996), are incorporated into the 
code. The t— u) model has both high and low Reynolds 
number formulations. The Cebeci-Chang(1978) rough- 
ness model is in both algebraic turbulence models. Re- 
sults are given for only one of the algebraic models, 
Chima et al. ( 1 993) . The effect of surface roughness was 
found to be relatively the same for both. 

C-type grids were used, and were generated using 
the method of Arnone et al.(1992). Grid parameters, 
especially the near wall spacing, were varied to ensure 
that the presented results were obtained with grids of 
sufficient resolution. A typical grid was 257 x 49 with a 
maximum near wall normalized grid spacing less than 
one. The turbulence models had differences in heat 
transfer sensitivity to near wall grid spacing. This is 
discussed in the comparison of results. The heat trans- 
fer coefficients are determined by the fluid temperature 
gradient at the wall 

Equivalent roughness height. For the algebraic and two- 
equation turbulence models roughness enters into the 
analysis by means of an equivalent roughness height, 
/*eq The equivalent roughness height accounts for vari- 
ations in the spatial distribution of roughness. Sev- 
eral different correlations have been proposed to obtain 
the h-E q from the roughness. These correlations require 
that roughness be characterized as having a geometric 
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Fig. 1 Comparisons of squlvalsnt height ratio con-stations. 

shape. Figure 1 compares different correlations for ob- 
taining the equivalent roughness height as a function of 
the roughness density parameter, A. These correlations 
were developed from data using deterministic, and not 
random, roughness. This figure is primarily illustrative. 
Direct comparisons cannot be made among the differ- 
ent correlations for roughness height because each uses 
a different definition for A. However, it does illustrate 
the large range of equivalent height ratios predicted by 
different correlations available in the literature. 

Six roughness traces were obtained from the vane 
surface tested by Boyle et. al.(2000). Figure 2 shows 
a typical roughness trace. From each trace the RMS 
height, /jrmSi w r as calculated The average value was 
111/zm, and the one sigma variation among the six 
traces was ±18//m. The profiles in the traces were then 
analyzed assuming the roughness was either cones or 
hemispheres. Table I gives the equivalent height ratios 
obtained using the different correlations for each of the 
six traces. If each roughness trace had a different value 
of /irms i but all were geometrically similar, there would 
be no variation in Aeq / ^rms- Results are shown for two 
model geometries; cones and hemispheres. Results are 
also shown for two assumptions as to what constitutes a 
peak in the roughness profile. In one definition, a peak 
is given as a change in the profile slope which occurs 
above the mean. In the other definition, a peak has to 
exceed the mean by Arms ■ The second assumption was 
made because it was felt that peaks in the shadow of 
higher peaks might not be as influential with respect 
to surface roughness effects. Compared to the first as- 
sumption, the second assumption leads to fewer, but 
taller peaks, spaced further apart. The results in Ta- 
ble I show that the second assumption results in only 
a small decrease of about 10 to 20% in the equivalent 
height ratio, Aeq/Arms- 



Fig. 2 Typical rougftrms trace. 

Even after accounting for different definitions of A, 
there is a significant variation for the equivalent height 
ratio among the correlations. When the roughness is 
modeled as cones, the Sigal and Dan berg( 1990) corre- 
lation gives an equivalent height ratio about half that of 
the other correlations. For conical models there is less 
variation among the Dvorak(1969), Dirling(1973), and 
W aigh and Kind(I998) correlations. For hemispheri- 
cal models the Dvorak, and Waigh and Kind correla- 
tions are significantly higher than the other two corre- 
lations. In this table there are two values of the equiv- 
alent height ratio for each model geometry and trace. 
The roughness traces were analyzed for skewness and 
kurtosis. The negative value of kurtosis indicates a 
bumpy rather than spiky profile, Dagnall(1986). This 
would favor hemispheres over cones as the model ge- 
ometry to determine the equivalent height. However, 
w 7 hen roughness elements are modeled as hemispheres, 
rather than cones, the slope of the roughness is not con- 
sidered. When modeled as hemispheres, the height and 
radius are equal. The average slope was calculated to be 
greater than 45 degrees. Also, hemispheres have greater 
frontal area than cones. For these reasons it is felt that 
conical elements were appropriate for determining the 
equivalent roughness height. The average equivalent 
height ratio, Arq/Arms, was 2.1 for the widely spaced 
conical roughness elements. 

Even though a single roughness height is used in 
the analysis, the normalized roughness height, A + varies 
around the vane in accordance with the local flow 7 con- 
ditions. 

h + = h Eq U FSP y/cJ/2/fi 

The freestream velocity, J7 fSi and p are obtained from 
pressure distribution using the isentropic relationships. 
When Cf/2 is zero, A + is zero, independent of the ac- 
tual roughness height. 
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Table I. Equivalent height ratio, AeqA 



Equivalent Height Roughness Model 

Trace 

/'rms mm 


!■ 

Dvorak 


Dirling 

Waigh A Kind 







Cone Model 





A 

ft eq Arms 

mm 

ft eq /Arms 

A 

ft eq Arms 

A 

IKZSQ3SHS 

1 

0.128 

45.1 

1.95 

11.5 

4.31 

10.0 

3.43 

27.1 

3.84 

2 

0.131 

44.6 

1.84 

10.9 

4.25 

9.9 

3.28 

27.3 

3.53 

3 

0.107 

47.1 

1.96 

12.1 

4.32 

10.3 

3.49 

28.2 

3.92 

4 

0.09T 

72.3 


18.8 

2.13 

12.7 

1.87 



5 

0.110 

87.2 


21.3 

1.79 

13.8 

1.57 

53.4 


6 

0,085 

152.4 


42.2 

0.90 

18.8 

0.95 

87.7 


Avg. 

0.111 


1,28 


2.95 


2.43 

■1 

2.67 






ft PEAK > 0 - 

Hemisphere 





A 


HI 

ft EQ/ ft RMS 

A 


IHBH 

ft eq Arms 

1 

0.128 

20.0 

■ 

6.6 

8.12 

8.1 

5.14 

4.9 

9.28 

2 

0.131 

18.1 


5.9 

8.46 

7.7 

5.25 

4.4 

7,15 

3 

0.107 

21.3 

5.56 



8.3 

5.16 

5.2 

11.05 

4 

0.097 

33.2 

2.51 

10.9 

3.93 


2.74 

19.5 

4.63 

5 

0.119 

35.4 

2.25 

11.7 

3.56 


2.51 

20.7 

4.22 

6 

0.085 

80.7 

0.84 

26.6 

1.52 

16.3 

1.26 

47.3 

1.99 

Avg. 

0.111 


3.80 

i 

5.61 


3.68 


6.39 





h 

PEAK > ft RMS 

- Cone Model 





A 


MM 


hh 


HSU 

ft eq Arms 

1 

0.128 

68.4 

| ■ 

17.5 




41.0 

3.14 

2 

0.131 

61.1 


15.0 


11.6 


37.4 

3.47 

3 

0.107 

53.8 


13.8 

4.79 


3.96 

32.1 

4.41 

4 

0.097 

98.6 


25.6 

1.79 

14.9 

1.67 

58.7 

1.79 

5 

0.119 

162.2 


39.7 

1.14 

18.8 

1.12 

99.2 


6 

0.085 

189.5 


52.5 

0.85 

21.0 

0.93 

109.0 

1.02 

Avg. 

0.111 


l.n 


2.66 ; 


2.31 


2.49 





h 

PEAK > ft RMS 

- Hemisphere 





A 

Aeq/Arms 

A 

Aeq/ ft RMS 

A 

ft eq Arms 

A 

ft eq Arms 

1 

0.128 

30.4 

4.10 

10.0 

6.32 

10.0 

4.33 

17.8 

7.37 

2 

0.131 

24.8 

5.37 

8.2 

8.01 

9.0 

5.28 

14.6 

9.11 

3 

0.107 

24.5 

6.02 

8.0 

8.94 

5.9 

5.86 

14.2 

10.15 

4 

0.097 

45.3 

2.05 

14.9 

3.31 

12.2 

2.44 

26.5 

4.04 

5 

0.119 

65.8 

1.29 

21.7 

2.27 

14.7 

1.80 

38.6 

2.89 

6 

0.085 

100.4 

0.76 

33.1 

1.43 

18.1 

1.23 

58.9 

1.91 

Avg. 

0.111 


3.27 


5.05 


3.49 


5.91 


Algeb raic turbulenc e m odel. The Cebeci-Chang(1978) 
roughness model increases the mixing length to aecoun t 
for roughness. The distance increment is given by: 

Ay+ =0.9(Vh+ - h + exp-° 167h+ ) 

and 

Ay=Ay + v/U FS py/cJ/2 

In the algebraic models the increment in y is only ap- 
plied in the inner region. The turbulent eddy viscosity 
in the inner region, fit, i is given by: 

IH,i = p(dU/dy) [«(y + Ay)(l - exp-^ ++ ^ A+ )f 


When fully turbulent flow is not assumed, the transi- 
tion model does not account for roughness. There is no 
mechanism in the analysis to allow for roughness effects 
when the flow is calculated to be laminar. 

k — u; turbulence model. In this model roughness influ- 
ences the solution through the boundary condition on 
u/. At the vane surface 

uo = max ^win, max(2500/(ft + ) 2 , 100 /h+)dU/dyj 

and 

WIN = *IN 

At the surface k 0 — 0. 
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Table II. Summary of test conditions, 



Ideal Axial chord Re No. xiO e 

M* 

Inlet 

Exit 

0.3 

0.0055 

0.041 


0.0086 

0.061 


0.0157 

0.101 


0.0278 

0.202 


0.0405 

0.303 

0.7 

0.0317 

0.236 


0.0521 

0.394 


0.0792 

0.590 

0.9 

0.0307 

0.240 


0.0555 

0.444 


0.0826 

0.665 



Fig. 3 Computational domain 

The inlet u ; decreases as the inlet turbulence intensity 
decreases, and as the length scale increases. All of the 
comparisons with data examined herein had a large in- 
let length scale when the inlet turbulence was high. At 
low Twin and large /in the wall boundary condition on 
uj is determined by h + . When win is small and A + > 25, 
u>o — 100 /h+dU/dy. When this condition applies flow 
is likely to be fully rough, where C//2 is independent 
of Reynolds number. Increasing Re at a constant Mach 
number causes h + to increase, which results in u>o de- 
creasing nearly in proportion to Re increasing. Also, 



Fig. 4 Predicted pressure coefficient. 

if the exit Mach number increases at constant Re, 
dU/dy increases, and ujq increases. Increasing vane exit 
Mach number at constant Re has the same effect on u?o 
as decreasing Re at constant exit Mach number. The 
difference between the high and low Reynolds number 
formulations of the k — w models is that the low formu- 
lation modifies the production and destruction terms by 
factors which are functions of the turbulent Reynolds 
number. Details of the implementation of the high 
and low Reynolds number formulations are gi\en by 
Chima(1996). For smooth surfaces, the low Reynolds 
number formulation, denoted as Lkw, is prefered. How- 
ever, it was found that this formulation predicted tran- 
sition further downstream than was seen in the data. 
The high Reynolds formulation, denoted as Hlw;, gave 
transition behavior more consistent with the data. 

DISCUSSION of RESULTS 
Stator data results 

The range of test conditions is given in Table II. Com- 
parisons are made with data for three or five Reynolds 
numbers at each of three exit Mach numbers. For each 
test condition shown in Table II data were obtained at 
two inlet turbulence levels. 

The vane, along with the computational domain 
is shown in figure 3. The grid shown in this figure is 
not the actual grid, since many grid lines have been 
removed to clarify the figure. The figure shows a vane 
with 80 degrees of turning. 

Figure 4 shows the pressure distribution around 
the vane at the three exit Mach numbers. The pressure 
coefficient, C p , is the difference between the inlet total 
pressure and the local static pressure divided by the 
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a) Comparison of Cebeci-Chang results with data. 



0 • — — - ‘ 1 ■ ' • 

-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

Surface distance s/C, 


b) Comparison of Hk« results with data. 



c) Comparison of Lk« results with data 




Surface distance, s/C, 
e) Comparison of Hkw results with data. 

Fig. 5 Nusseit number comparisons at M } =0.3, no turbulence grid. 

exit ideal dynamic pressure. The pressure distribution 
at each Mach number is for data obtained at a midrange 
Reynolds number. The variation in C p with Reynolds 
number at constant M o was found to be small, as was 
the variation in C p with surface roughness. On the 
suction surface C P shows a rapid increase; the rate of 
which decreases as the Mach number increases. At the 
maximum C p the static pressure is less than the vane 
row exit pressure. The adverse pressure gradient which 
follows assures that, even at low turbulence intensity, 
for most test cases the flow over the remainder of the 
suction surface will be turbulent. The adverse pressure 
gradient increases as the Mach number increases. The 
pressure surface has low 7 velocities over the first 40% of 
the surface followed by strongly accelerating flow 7 . The 
same analysis as is used here for the heat transfer pre- 
dictions w T as used to compare predicted and experimen- 
tal pressures for a smooth vane of the same geometry. 
Boyle et al.(1998) show r ed good agreement between pre- 
dicted and measured C p for the smooth vane. 

He at tr ansfer comparis ons. Heat transfer is given in 
terms of Nusseit number, using C x as the reference di- 
mension, and the thermal conductivity at the inlet total 
temperature. Calculations w 7 ere done with the same av- 
erage gas-to-wall temperature ratio as for the data. Be- 
cause the temperature differences were small, the Nus- 
seit number was based on the difference between the 
surface and the local adiabatic recovery temperature. 
The local recovery temperature w 7 as determined from 
the local static pressure, the iseri tropic relationships, 
and a recovery factor of 0.9. 

For each test condition comparisons are show r n using 
one of the three models. Because the roughness height 
is not known precisely, and to show 7 the sensitivity of 
predicted heat transfer to variations in roughness, pre- 
dictions are showrn for a range of roughness heights. 
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Under some conditions the sensitivity was high, while 
for other conditions the sensitivity was small. Predic- 
tions are shown in order of increasing Mach number. At 
each Mach number, heat t ransfer comparisons are made 
at various Reynolds numbers for low turbulence inten- 
sity, and are followed by predictions at the high tur- 
bulence intensity. Predictions are shown for the entire 
vane surface, not just those regions where heat transfer 
data were available. While eac h test condition shows 
predictions for only a single model, successive compar- 
isons use different models. In this way, the effects of 
different models can be compared for similar, though 
not identical test conditions. The three models used 
are identified as: (1) the Cebeci-Chang model, (2) the 
high Reynolds number k — model, (Hku/), and (3) the 
low Reynolds number k — u model, (Lku;). 

Figure 5a shows comparisons for A/ 2 — 0.3, J?e 2 = 
0.041 x 10 6 , and no grid using the Cebeci-Chang rough- 
ness model. The heat transfer at this test condition 
is most likely to resemble smooth vane heat transfer. 
At sjC x near 2.6 on the suction surface, transition is 
seen in the data. Prior to this location, the transi- 
tioning analysis, which in reality is for laminar flow, is 
somewhat higher than the data. The transition model 
did not account for roughness, so the flow remained 
laminar. The predicted suction surface heat transfer 
nearly doubled when assuming fully turbulent flow from 
the leading edge. Even for an equivalent roughness 
height twice the expected equivalent height, the surface 
roughness is not large. This is expected, since at this 
Reynolds number the ratio of /* + //?eq is a minimum. 

The experimental data, which were obtained from an 
infrared camera, are shown as up and down triangles. 
In the forward part of the vane, data from each camera 
view are indicated by different symbols. Moving from 
the leading edge at s/C x = 0.0 to the pressure surface 
both the analysis and the data show a very rapid de- 
crease in Nusselt number. The experimental data are 
shifted further along the pressure surface of the vane. 
This may be partially due to inaccuracies in mapping 
the camera view onto the vane surface. No adjustment 
was made in the mapped coordinates to achieve better 
positional agreement in the leading edge region. While 
this could be justified, the offsets seen in this and sub- 
sequent figures allow clearer comparisons between the 
data and analyses. In the leading edge region the ex- 
perimental heat transfer is higher than predicted. The 
predicted stagnation point Frossling number is less than 
one, while, as expected the experimental value is close 
to one. In the analysis h + varies in proportion to yf. 
Near the leading edge exceeded the suction surface 
value at s/C x = 2.6, where the data show transition. 


Figure 5b compares experimental Nusselt numbers 
with the Hku; predictions for the next highest Reynolds. 
Experimentally, suction surface transition mowd for- 
ward with increased Reynolds number. This turbu- 
lence model shows smooth surface transition closer to 
the leading edge than the data. Increasing roughness, 
and thus the boundary value for w, moved the predicted 
transition location forward. Predicted Nusselt numbers 
are more sensitive to roughness variations than for the 
previous case. This is expected because /?+ increases 
with Reynolds number for a given roughness height. 

Data and predictions are shown for the Lku; tur- 
bulence model for i?e 2 = 0.101 x 10 6 in figure 5c. Pre- 
dictions are shown for inlet turbulence intensities of 3 
and 5%. The figure shows that transition occurred only 
for an inlet turbulence of 5%, and Z^eq/Cx = 0.01. No 
predictions are shown for the experimental inlet turbu- 
lence of 1%. Steady state solutions were not obtained at 
experimental inlet turbulence levels due to oscillations 
caused by vortex shedding. There is poor agreement 
between the analysis and the data. In the leading edge 
region the analysis underpredicts the Nusselt numbers. 
On the suction surface the data show transition begin- 
ning at s/C x close to 1.0. The post transition Nu of 
about 300 is consistent for turbulent flow with the max- 
imum level of 200 seen in figure 5b, for a 50% increase 
in Reynolds number. 

Figure 5d shows a second Nusselt number compar- 
ison using the Cebeci-Chang roughness model, assum- 
ing fully turbulent flow. Since the local values of h + 
are greater at the higher Reynolds number, the rela- 
tive effects of surface roughness are also greater. The 
Cebeci-Chang model predicts the heat transfer in the 
leading edge region well, but underpredicts the effects 
of surface roughness on the rear portion of the suction 
surface. As in part (a) the fully turbulent analysis gives 
very high suction surface heat transfer near s/C x = 0.5. 
The data show rapidly increasing heat transfer prior to 
this location. However, the peak experimental value 
could not be determined, because the infrared camera 
could not see further along the suction surface. 

The Hku; predictions are compared with data for 
the highest Reynolds number case at M 2 = 0.3 in figure 
5e. The predicted effect of surface roughness on both 
the pressure and suction surfaces is large. The analysis 
agrees well with the data in the leading edge region, 
and for the rear of the suction surface. On the suction 
surface the prediction shows a rapid transition occur- 
ring close to the leading edge for the rough surface. The 
data indicate, however, that transition is not complete 
until s/C x = 2. Even for the smooth surface, the analy- 
sis shows transition occurring closer to the leading edge 
than does the data. 
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a) Comparison of Lk«> results wfth data. 



b) Comparison of Cebeci-Chang results with data. 
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c) Comparison of Hk« results with data. 



d) Comparison of Lka> results wfth data. 



Surface distance, a/C, 

e) Comparison of Cebeci -Chang resufts with data. 

Fig. 6 Nussett number comparisons at M ? =0.3, turbulence grid. 

Next, comparison with data are shown for the same 
exit Mach number, but for a higher turbulence level 
caused by the presence of the grid. 

Figure 6a shows results for the lowest Reynolds 
number using the Lku; model. The influence of surface 
roughness on the predicted start of suction surface tran- 
sition, and on the post transition Nusselt number level 
is clearly seen. The data show a two step transition 
process. However, Boyle et al.(2000) pointed out that 
at the lowest Reynolds number there is the possibility 
that, very close to the trailing edge, the Nusselt num- 
ber is in error because heat loss across the thickness of 
the vane. This source of possible error decreases as the 
Reynolds number increases. 

The Cebeci-Chang model predictions are compared 
with data in figure 6b for the second lowest Reynolds 
number. In the leading edge region there is a significant 
difference in heat transfer between the two views. If the 
average of the two values from the two views is taken, 
the agreement in the overlap region is reasonably good. 
Since the estimated equivalent height ratio, Aeq /C r is 
0.005, the analysis agrees well with the experimental 
data for most of the reaT portion of the suction surface. 
For the part of the pressure surface for which data are 
available, the analysis correctly predicts the rapid de- 
crease in heat transfer moving away from the stagnation 
point. 

Comparison with data for the middle Reynolds 
number are shown for the HkuJ model in figure 6c. The 
predictions are lower than the experimental data for 
the rear portion of the suction surface. At the leading 
edge the heat transfer level is well predicted based on 
the suction side data. Much of the suction side data 
show r decreased heat transfer from the leading edge. A 
minimum is reached near s/C x — 0.15. Then the heat 
transfer increases rapidly. None of the analyses show 
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a) Comparison of Cebeci-Chang results with data. 



b) Comparison of Hkw results with data. 



a) Comparison of Lkw results with data. 
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b) Comparison of Cebeci-Chang results with data. 
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c) Comparison of Lk<o results with data 
Fig. 7 Nusseit number comparisons, M 2 =0.7, no turbulence grid. 

this suction surface trend. In the analyses the Nus- 
seit number either remains fairly constant, or, at higher 
Reynolds numbers increases very rapidly after the stag- 
nation point. 

The Lku; predictions shown in figure 6d are signif- 
icantly lower than the data for the rear portion of the 
suction surface. The experimental heat transfer in the 
leading edge region is higher than the prediction. 


c) Comparison of Hko results with data 
Fig. 8 Nusseit number comparisons for M 2 =0.7, turbulence grid. 

Comparisons using the Cebeci-Chang model are 
shown in figure 6e for the highest Reynolds number test 
at M 2 = 0.3. In the leading edge region the agreement 
of the fully turbulent prediction at /?eq/Ct = 0.005 
with data is good. The Nusseit numbers at the rear of 
the suction surface are significantly underpredicted. On 
the suction surface the Cebeci-Chang roughness model 
shows a small decrease in Nusseit numbers when the 
roughness is increased from 0.005 to 0.01. This is 
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probably due to the increment in the mixing layer from 
roughness being only applied in the inner boundary 
layer region. For very high roughness, (in terms of A + ), 
this model may not give the appropriate increase in heat 
transfer with increasing roughness. This may also be 
true of the decrease in blade row efficiency with rough- 
ness. In terms of the predictions, the largest effect, of 
roughness was to cause fully turbulent flow. The in- 
creased heat transfer on both suction and pressure sur- 
faces due to assuming fully turbulent flow was greater 
than the increase due to surface roughness. 

Raising the exit Mach number from 0.3 to 0.7 caused 
a large increase in the minimum exit Reynolds number. 
In addition to more than doubling the exit velocity, the 
minimum inlet total pressure increased. Consequently, 
the minimum Re 2 at Mo = 0.7 was nearly six times 
greater than the minimum Re 2 at Mo — 0.3, and was 
80% of the maxim urn Ren at, Mo = 0.3. 

Figure 7a compares experimental and predicted Nus- 
selt numbers for the lowest Reynolds number tested at 
A/ 2 — 0.7. These results are for the no grid case, and 
show the Cebeci-Chang predictions. If Mach number 
did not affect the heat transfer, it is expected that the 
experimental data would be nearly the average of that 
shown in figures 5d and 5e. In the leading edge region, 
where the local Mach number is low, this is seen to be 
true. However, on the rear portion of the suction sur- 
face, where the local Mach number is close to the exit 
value, the variation of Nusselt number with surface dis- 
tance is different. The higher Mach number data in 
figure 7a show a more negative slope over the rear por- 
tion of the suction surface. Figure 7a shows that the 
analysis agrees well with the data in the leading edge re- 
gion. However, it is much lower than the measurements 
for the rear of the suction surface. 

Figures 7b and 7c show comparisons with data 
for the two versions of the k — w model. In the lead- 
ing edge region, both models show satisfactory agree- 
ment. On the rear of the suction surface both models 
are in relatively good agreement with the data. The 
Hk u model appears to have somewhat better agree- 
ment. However, both models agree with the data to 
the same degree. Either turbulence model, when com- 
pared with the other Reynolds number data, showed 
the same degree of agreement with the data. The re- 
sults in these two figures also illustrate that, continually 
increasing the roughness height results in only small 
changes in the predicted heat transfer. 

Comparing the results in figures 7 and 8 show 
the effect of an increased turbulence level. There is 
no consistent experimental heat, transfer variation with 
turbulence intensity for the rear of the suction surface. 



a) Comparison of Cebeci-Chang results with data. 



b) Comparison of Hkco results with data. 



-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

Surface distance, s/C„ 

c) Comparison of Lko) results with data. 

Fig. 9 Nusselt number comparisons, M,=0.9, no turbulence grid. 

But, the variation is small; within the measurement un- 
certainty at each Reynolds number. In the leading edge 
region the variation is more consistent, and shows a heat 
transfer increase with higher turbulence intensity. 

Comparing predictions to each other for the 
cases presented in figures 7 and 8 illuminate differ- 
ences among the turbulence models. The Cebeci-Chang 
model shows a smaller increase in heat transfer with 
roughness than do either of the two k — u? models. The 
Hku> model is in good agreement with the measured 
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a) Comparison of Lkw results with data. 



heat transfer for the rear of the suction surface. Figure 
8, and to a lesser extent figure 7, show that for the 
forward portion of the vane, on either side of the leading 
edge, all models agree reasonably well with the data. 
This occurs when, in the overlap region, the average 
from the two views is used. 

The results shown in figure 8c illustrate a difficulty 
that was observed with the Hku> model. The predictions 
show a maximum Nusselt number at an intermediate 
roughness height. This was surprising, and the cause 
for this was found to be sensitivity of the HIw model 
to the near- wall grid line spacing. 

Lastly, comparisons are shown in figures 9 and 10 
for Mo — 0.9. The percentage change in exit Reynolds 
number is not as great as the the percentage change 
in exit Mach number. At high Mach numbers density 
decreases with increasing Mach number, resulting in 
smaller Reynolds number changes. Comparing results 
in figures 7 and 9 show the effect of a Mach number 
variation for no grid-low turbulence cases. Comparisons 
between figures 8 and 10 show the Mach number effects 
for high turbulence. 



c) Comparison of Hkco results with data. 



Surface distance, s/C x 


d) Comparison of Lko> results wfth data. 

Fig. 10 Nusselt number comparisons for M 2 =0.9, turbulence grid. 

The smooth surface results for fully turbulent flow 7 
in figure 9a show 7 a sharp decrease in suction surface 
Nusselt number beginning at s/C x = 1. While this 
was seen to some extent for the transitioning predic- 
tion in figure 7a, it was not seen for the comparable 
fully turbulent prediction. Figure 4 shows that the ad- 
verse suction surface pressure gradient is steeper for the 
higher exit Mach number. The steeper gradient leads to 
a more rapid thickening of the boundary layer, leading 
to the decrease in heat transfer. Although not shown in 
the figures, the two k — u model results showed similar 
Nusselt number behavior to that seen in figure 9a for 
the algebraic model. The magnitude of the dip in heat 
transfer, however, was less for both k — u> models. The 
dip in heat transfer is less when the analyses are for a 
rough surface. The experimental data in figure 9 also 
exhibit this behavior, but so do the data in figure 7. 


NASA/TM — 2000-210219 


11 



Table III. Summary of rotor test conditions. 



Ideal Axial chore 

Re No. x 10 -e 

An 

Inlet 

Exit 

50 

0.37 

0.56 


0.28 

0.42 

36 

0.23 

0.42 


0.12 

0.23 


In the leading edge region figure 9 shows that, for all 
three Reynolds numbers, the analysis agrees reasonably 
well with the data. There appears to be little reason 
to prefer one turbulence model over the other based on 
the data in this figure alone. 

The comparisons in figure 10 for Mo — 0.9 show 
results very similar to those in figure 8 for Mo = 0.7. It 
was found that, of the two k — w turbulence models, the 
model without a low Reynolds number effect (the Hkw 
model), was more sensitive to near wall spacing than 
was the Lkw model. Even though the grid had a near 
wall spacing of yf = 1 or less, the Hkw model results 
changed when the near wall spacing w T as reduced. The 
HKw predictions for the highest Reynolds number case 
shown in figure 10c became almost identical to those 
shown in figure lOd for the Lkw model when the near 
wall grid spacing was reduced. The Lkw model and the 
Cebeci-Chang model heat transfer predictions did not 
change with reduced near wall spacing. 

The term win used as part of the criteria for deter- 
mining wq is not part of the standard k — w turbulence 
model, Wilco x(1994). It is most lilely to have an in- 
fluence at high Reynolds numbers and high turbulence 
intensities. For the cases examined in this study, omit- 
ting win from the criteria for determining wo did not 
affect the heat transfer predictions. 

Rotor blade results 

BIair(1994) presented both full span and midspan 
rotor heat transfer distributions at different Reynolds 
numbers and incident flow angles. Measurements were 
made for both a rough and smooth blade. In contrast 
to the stator data, these rotor blade data had mea- 
surement over the entire surface. Table III describes 
the test conditions for which comp arison s are made. 
For the rough surface a screened grit of 0. 055mm was 
applied to the surface. This grit size is 0.0041 when 
normalized by the reference chord, C x Additional in- 
formation regarding the roughness height was not avail- 
able. Roughness calculations were done using multiples 
of this height . A value of /*eq /C x = 0.0041 corresponds 
to an equivalent height- to- RMS value, Aeq/Arms, of 
approximately two. 



Surfac* distance, s/C I 


Fig. 11 Predicted pressure coefficient for rotor of Blair(1994). 

The pressure coefficient distributions are shown in 
figure 11 for the two inlet flow angles. At the higher in- 
let flow angle of 50 degrees there is an overspeed close 
to the leading edge of the suction surface. This causes 
rapid transition. At this inlet angle the suction surface 
will be almost completely turbulent . At the lower inlet 
flow angle of 36 degrees, there is less flow turning. Con- 
sequently, the pressure difference across the rotor is less. 
At this inlet angle the flow accelerates more uniformly. 
The smooth surface boundary layer remains laminar 
until the peak pressure coefficient, Cp } is reached at 
s/C x of approximately 0.9. 

At an inlet angle of 50 degrees, the pressure surface 
velocities are very low for nearly the first half of the sur- 
face distance. This is followed by a strong acceleration. 
The smooth surface pressure surface boundary is likely 
to be laminar, especially at the lower Reynolds num- 
ber. For an inlet angle of 36 degrees, there is a pressure 
surface overspeed. Unless the flow relamin arizes, the 
pressure surface boundary layer will be turbulent. 

Figure 12 compares measured and predicted Nus- 
selt numbers for the rotor at the four test conditions. 
In each case two rough surface predictions are shown. 
They are for roughness height ratios. Aeq/C* of 0.004 
and 0.008. The two values shown illustrate the sen- 
sitivity of the predictions to roughness height. The 
higher value represents an estimated upper bound on 
the roughness height. 

The comparisons in figure 12a are for An = 50°, 
and an exit Reynolds number of 0.56 x 10 6 using the 
Hibw model. The agreement with both the rough and 
smooth surface data is good for almost the entire rotor 
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a) Comparison of Hkw results withdaia. 



b) Companion of LM«> results with data. 
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c) Comparison of Cebeci-Chang results with data. 



Fig. 12 Comparisons with data of Blair(1994) . 


surfaces. Only for the first third of the pressure sur- 
face does the analysis underpredict the rough surface 
heat transfer data. Here again, the model shows 
decreased suction surface heat transfer with increased 
roughness. 

The predictions shown in figure 12b are for the L ku> 
model. The results are for the lower Reynolds number , 
but still at the design inlet flow angle. The agreement 
with data is good for both the smooth and rough sur- 
faces. The predictions are lower than the rough surface 
data for the forward portion of the pressure surface, 
and high than the data midway along the suction sur- 
face. The Lku> model shows increased heat transfer as 
the surface roughness roughness is doubled. This is in 
contrast to the H model results shown in figure 12a. 

The analysis does not agree as well with the 
data for the two off-design inlet flow angle cases. The 
Cebeci-Chang model results shown in figure 12c under- 
predict the effect of roughness. At a roughness height, 
/ieq fC x = 0,004 the Cebeci-Chang model shows a sig- 
nificant increase due to roughness in predicted suction 
surface heat transfer. However, the model predicts al- 
most the same heat transfer as for a smooth surface, 
when the roughness height is doubled. This suprising 
result is consistent with the heat transfer results ob- 
tained using this model for the stator vane. At high 
values, this model showed a peak in surface heat trans- 
fer, as the roughness height increased. The analysis is 
in reasonably good agreement with the smooth surface 
data for the suction surface. However, it underpredicts 
the pressure surface heat transfer close to the leading 
edge, and overpredicts the heat transfer closer to the 
trailing edge. The smooth surface transitioning model 
gives higher heat transfer over much of the pressure sur- 
face, than does the fully turbulent analysis. This occurs 
because the transitional model augments the eddy vis- 
cosity in the laminar region to account for freest ream 
turbulence effects. The fully turbulent model does not 
augment the eddy viscosity to account for freestream 
turbulence effects. 

The Lkuj model results shown in figure 1 2d for 
the lower Reynolds number show that for the suction 
surface this model accurately predicts heat transfer for 
the rough surface. Only for the rear half of the pressure 
surface does this model agree with the experimental 
rough surface data. The results for the forward half of 
the pressure surface are similar to those shown in figure 
12a for the design inlet flow angle. However, the degree 
of agreement with data is poorer for the lower inlet flow 
angle. For the smooth surface the L kw model predicts 
suction surface transition before it is seen in the data. 
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Figure 12c shows that the transition model used with 
the algebraic turbulence model gives better agreement, 
with the smooth blade data. 

Overall, the results for the rotor blade comparisons 
show the same trends as for the stator vane. Where 
roughness causes transition dose to the leading edge, 
the H km model gave best agreement with data. If 
the momentum thickness Reynolds numbers were high 
enough to cause smooth surface transition close to the 
leading edge, the L km model gave good agreement for 
roughness effects. The Cebeci-Chang roughness model, 
while it does show increased heat transfer for rough sur- 
faces, underpredicts the effects of surface roughness. 

CONCLUSIONS 

The primary conclusion of this work is that the high 
Reynolds number formulatiai of the k — m turbulence 
model results in the best agreement with the experi- 
mental data. It is best in the sense that it provides a 
conservative estimate for the effect of roughness on heat 
transfer. The high Reynolds number formulation gave 
early transition at moderate Reynolds numbers and low 
turbulence intensities. This was consistent with the ex- 
perimental results. The low Reynolds number formula- 
tion showed transition like behavior at low turbulence 
intensities, which was inconsistent with the data. 

Predictions made using the Cebeci-Chang turbulence 
model showed a heat transfer increase with roughness. 
When this model w as used for a rough surface, the flow 
was assumed to be fully turbulent. This approach re- 
sulted in better agreement with the data than the low 
Reynolds number formulaticn of the k — m turbulence 
model. The heat transfer increase due to roughness was 
less with the Cebeci-Chang model than with the high 
Reynolds number k — m model. Comparisons with the 
data showed the high Reynolds number k — m predic- 
tions agreed better than the Cebeci-Chang model at 
higher Reynolds numbers. This is consistent with the 
observations of Boyle and Civinskas(1991). They re- 
ported lower than measured heat transfer for some test 
cases using the Cebeci-Chang roughness model.” 

The choice of model for rough surface heat transfer 
predictions is strongly influenced by accurate knowl- 
edge of the equivalent roughness height. In addition 
to the variation in roughness height at various sur- 
face locations, there was a large variation in equivalent 
height among various models for predicting the equiva- 
lent height . Reducing the variation among correlations 
for equivalent height is as important as improving the 
turbulence model for rough surface heat transfer pre- 
dictions. 
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